Squeezing and entanglement of matter-wave gap solitons 
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We study quantum squeezing and entanglement of gap solitons in a Bose-Einstein condensate 
loaded into a one- dimensional optical lattice. By employing a linearized quantum theory we find 
that quantum noise squeezing of gap solitons, produced during their evolution, is enhanced compared 
with the atomic solitons in a lattice-free case due to intra-soliton structure of quantum correlations 
induced by the Bragg scattering in the periodic potential. We also show that nonlinear interaction 
of gap solitons in dynamically stable bound states can produce strong soliton entanglement. 
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In the past two decades, the active research into quan- 
tum properties of nonlinear many-particle systems has 
enabled control and engineering of quantum noise in non- 
linear optics P]. Quantum noise reduction of optical 
signals below the shot-noise level, referred to as uncer- 
tainty squeezing, and quantum entanglement of inter- 
acting optical pulses are at the core of the applications 
of non-classical light in quantum interferometry, preci- 
sion measurement, and quantum information process- 
ing. Non-spreading optical collective excitations, soli- 
tons, supported by dispersive nonlinear media, proved 
the best candidates for experiments on quantum noise 
reduction and entanglement due to their robust dynam- 
ics and scattering properties [J . 

It has been recognized that general methods of 
quantum noise squeezing and entanglement developed in 
quantum optics could apply to other nonlinear bosonic 
fields, such as weakly interacting ultracold atoms in a 
Bose-Einstein condensate (BEC). Consequently, a num- 
ber of theoretical proposals wereput forward lEHill] and 
experiments carried out 0, 0, on generating macro- 
scopic entangled number-squeezed states in BEC. Pro- 
duction of quantum correlations in the macroscopic quan- 
tum states relies on interactions between the atoms, the 
mechanism that is analogous to optical Kerr effect. 

Recently, gap solitons in a repulsive condensate, sup- 
ported by periodic potentials of optical lattices, have 
attracted a great deal of attention due to their con- 
trollable interaction and robust evolution uninhibited 
by collapse [1 IHIll IH 111. The existence of gap 
solitons is a unique property of nonlinear periodic sys- 
tems. Atomic gap solitons form due to combination of 
an inherent nonlinearity of BEC and the bandgap struc- 
ture of the matter-wave spectrum that can modify and 
ultimately reverse dispersion properties of the atomic 
wavepackets The current techniques for gap soli- 

tons generation suffer greatly from "technical" noise |0 . 
Provided this problem can eventually be overcome, gap 
solitons may represent an attractive high-density source 
for atomic interferometry, quantum measurements and 
quantum information processing with ultracold atoms. 



Understanding and control of quantum noise associated 
with atomic gap solitons is therefore a fundamental is- 
sue which so far has not been explored. However, the 
studies of quantum polariton solitons in a frequency dis- 
persive medium |15|| and the amplitude-squeezed opti- 
cal solitons in shallow periodic media described by the 
couple-mode model suggest that both strong disper- 
sion and bandgap spectrum of a nonlinear soliton-bearing 
system with periodically varying parameters can dramat- 
ically modify the properties of the soliton squeezing. 

In this Letter, we study quantum fluctuations of 
matter-wave gap solitons in an optical lattice and inves- 
tigate the band-gap effect on the quantum noise squeez- 
ing. We employ the soliton perturbation approach [T^ 
to analyze quantum fluctuations around the soliton so- 
lutions of the Gross-Pitaevskii (GP) equation with a pe- 
riodic potential. Using this approach, we demonstrate 
enhancement of quantum noise squeezing effect induced 
by the gap soliton evolution. Furthermore, we show that 
the existence of dynamically stable bound states of gap 
solitons provides a favorable environment for generating 
entangled soHton pairs. 

We consider an elongated cigar-shape BEC loaded into 
a one-dimensional optical lattice and described by the GP 
equation for the macroscopic wave function [T]| , 

a* 1 

= "2 9^ + ^("^^^ + .9id|*|'*, (1) 

where ^' is the normalized macroscopic wave function of 
the condensate, guj is the nonlinear interaction coeffi- 
cient, and V{x) — Vosin^(x) is the one-dimensional peri- 
odic potential. Stationary states of the condensate can be 
presented in the form, ^/{t, x) = ip{x) exp(— i/it), where /i 
is the chemical potential. In the linear, non-interacting 
limit, gi£) — > 0, the spectrum of matter waves has the 
characteristic band-gap structure 11]. The nonlinear lo- 
calization of matter waves in the form of gap solitons 
occurs in the gaps of the linear spectrum, and the fam- 
ily of the lowest-order gap solitons in the first (finite) 
spectral gap is shown in Fig. ^ as the dependence of the 
soliton norm, P = J ip'^dx, vs. chemical potential. 
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FIG. 1: (Color online) Left: Band-gap diagram for Bloch 
waves (bands are shaded), and the family of the gap solitons 
in the first finite gap for Vo ~ 4.0. Right: Profiles of gap 
solitons with difi'erent chemical potentials; /i — 1.91, 3.0, and 
3.85, corresponding to the points a, b, and c, respectively. 

Since the degree of localization of gap solitons varies 
across the gap [see Fig.^a-c)], near the bottom edge of 
the gap, /i « /iQ, the weakly localized soliton profile is 
well described by the "envelope" approximation H^l, 

^P{x;fi)^AF{x)'i>{x;fio), (2) 

where <i>(x;/io) is the periodic Bloch state at the corre- 
sponding band edge, and F{x) is a slowly varying func- 
tion. The inset in Fig. [3 shows the oscillating wavefunc- 
tion ipix) of a gap soliton near the band edge, at = 1.91, 
together with the corresponding Bloch-wave envelope 
F{x). The envelope function F{x) is the solution of the 
lattice-free GP equation or nonlinear Schrodinger (NLS) 
equation with the effective anomalous diffraction and in- 
teraction energy modified by the lattice j2| • 

To study the quantum fluctuations of the gap soli- 
tons, we replace the 'classical' mean field described by 
Eq. ^ by the bosonic field operator ^P. Then, we 
use the linearization approach with the perturbed quan- 
tum field operator ip around the mean-field solution 'i>Q. 
By using the back-propagation method 19], we calculate 
the optimal squeezing ratio for gap solitons as a func- 
tion of the evolution time. The optimal squeezing ratio 
of the quadrature field for gap solitons, as a result of 
the homodyne detection scheme, can be defined as [T^ : 
R{t) = min[var(^'L(i)|'(/;(t))/var(*L(i)|V'(0))], where the 
inner product between the operator field, and the local 
oscillator profile, is defined as: 




where we use the mean-field soliton with an adjustable 
phase shift as the local oscillator profile, = 5*06*^ /P. 
By varying the phase of the homodyne detection, 0, we 
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FIG. 2: (Color online) Time evolution of the optimal squeez- 
ing ratio, R, for the gap soliton (solid), its near-band-edge 
envelope (dot-dashed), and the NLS soliton (dashed). Inset: 
profiles of the gap soliton found from Eq. (0 at = 1.91 and 
its envelope used in the approximation 

can determine the minimun value of the squeezing ratio 
of the quadrature component of the quantum field. When 
= 0, the in-phase quadrature component is detected, 
and when 6 = tt/2, the out-of-phase quadrature compo- 
nent is detected. The result of the calculations of R{t) 
for the gap soliton near the bottom edge of the spectral 
gap is presented in Fig.|21by a solid curve. 

The optimal squeezing ratio for gap solitons can be 
compared with that for a conventional NLS soliton co- 
inciding with the soliton envelope F{x) and evolving in 
a lattice-free anomalous diffraction regime (dashed line 
in Fig. [2J|. The NLS soliton of the same envelope is 
only weakly quadrature squeezed during propagation due 
to the nonlinearity of the matter-wave, whereas the gap 
soliton shows enhanced quantum noise squeezing. In con- 
trast, the near-band-edge envelope approximation of Eq. 
^ (dot-dashed line in Fig. |2Jl provides a good estimate 
for soliton squeezing in the initial stage of evolution. 

When the chemical potential of the gap soliton moves 
away from the band edge, the envelope approximation (0) 
becomes invalid. In Fig.lSfa), we show the dependence of 
the optimal squeezing ratio of a gap soliton at different 
values of the chemical potential inside the gap for differ- 
ent times. Towards the middle of the gap, the squeezing 
ratio improves with the chemical potential as the peak 
density of the gap soliton and hence the nonlinearity are 
increasing. Near the band edge the quantum noise of 
a gap soliton is squeezed in the quadrature field compo- 
nent, Fi g.EI b). Similar to the amplitude-squeezed Bragg 
solitons llal, the maximum squeezing angle approaches 
in-phase quadrature — > 0, when the chemical poten- 
tials of gap solitons move deeper inside the gap. However, 
number-squeezing is never achieved because we don't use 
an optimal profile for the local oscillator. For a fixed 
time and varying chemical potential, the best quadrature 
squeezing occurs around the middle of the gap. 

To understand the band-gap effect on the quantum 
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where Afij is the atom-number fluctuation in the j-th 
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FIG. 3: (Color online) (a) Optimal squeezing ratio, R, for the 
gap solitons with different chemical potentials, n, at different 
time, (b) Time evolution of the phase in the homodyne de- 
tection for the optimal squeezing ratios at different values of 
the chemical potential within the gap. 
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FIG. 4: (Color online) Quantum correlation spectra in the 
spatial (a;)-domain for gap solitons at different points within 
the gap and at the values of 9 corresponding to optimal 
quadrature squeezing: (a) 6 = 41.4°, (b) 9 = 9°, (c) 9 = 5.4°, 
and (d) 9 = 5.4°. Insets show the corresponding soliton com- 
ponents in the a;-domain. 



fluctuation of matter-wave gap solitons, we analyze the 
quadrature correlations between different spectral com- 
ponents of the gap soliton induced by its nonlinear evo- 
lution. The intra-soliton correlation coefflcients, Cij, are 
found by calculating the normally-ordered covariance, 



(3) 



slot Asj in the spatial (s 
domain: 



x) or momentum (s — k) 



An," 



ds[«'L(t,s)V't(t,s) + *2(t,s)V'(i,s)], 



The corresponding correlation spectra in the spatial (x)- 
domain for gap solitons at different points within the 
gap and for values of 9 corresponding to the maximum 
squeezing are shown in Figs, ^a-d), respectively. The 
near-band-edge gap soliton has a discrete correlation pat- 
tern "shaped" by the periodicity of the lattice, with 
the spatial components being weakly anti-correlated or 
uncorrelated [Fig. ^a)]. The anti-correlation of the 
gap soliton components explains its enhanced quadrature 
squeezing compared to the NLS soliton without a lattice. 

In the middle of the gap, the soliton is strongly spa- 
tially localized, and its spatial components are strongly 
anti-correlated [Fig. EIb,c)], which enhances the quadra- 
ture noise squeezing. However, near the top edge of 
the gap, the locahzation of the gap soliton degrades due 
to resonance with the Bloch state at the corresponding 
edge, and strongly correlated, noisy components asso- 
ciated with the periodic Bloch wave structure start to 
dominate in the correlation pattern [Fig. EJd)]. The 
squeezing is reduced near that edge of the spectral gap. 

The possibility of number squeezing for solitons can 
be assessed by analyzing the intra-soliton number corre- 
lations ed 6 = 0. In the momentum domain NLS soliton 
displays the well-known symmetric correlation pattern 
with the noisy, strongly correlated outer regions of the 
soliton spectrum jQ. For this reason, an efficient num- 
ber squeezing of the NLS solitons can be produced by 
spectral filtering that removes the noisy spectral com- 
ponents |2lj |. On the contrary, the number correlation 
spectrum of a near-band-edge gap soliton displays a pe- 
riodic pattern with regions of strong correlations at the 
Bragg condition, ki = kj — ±(2to + 1) (m is an integer) 
and uncorrelated off-diagonal regions. This correlation 
pattern persists through the entire gap until the noisy 
spectral components due to resonance with a Bloch state 
appear near the opposite band edge. Therefore the en- 
hancement of number squeezing of the BEG gap soliton 
with additional spectral filtering seems unfeasible due to 
spectral selectivity of the Bragg scattering. 

The possibility to achieve squeezing of atomic solitons 
paves the way to creating entangled macroscopic coher- 
ent states in BEG interference experiments. However, 
optical lattices also offer an alternative possibility for 
entanglement production. It is known that the optical 
lattice supports stationary states in the form of bound 
soliton pairs [Tlj . The nonlinear interaction between 
solitons in a pair could induce entanglement and hence 
nonseparability of the bound state. This entanglement 
can potentially be exploited in BEG atom number de- 
tection by implementing analogs of quantum nondemo- 
lition measurements based on entangled pairs of optical 
solitons l^^. Potential advantage of gap solitons com- 
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FIG. 5: (Color online) (a) Optimal squeezing ratio for in- 
phase (inset, solid) and out-of-phase (inset, dashed) bound 
soliton states (Vb = 4.0, \jl = 2.5); the two curves overlap, (b) 
Atom number correlation C12 for in-phase (solid) and out-of- 
phase (dashed) soliton pairs. 

pared to the lattice-free atomic solitons is the long-lived 
nature of the bound states which are pinned by the lat- 
tice and do not break-up as a result of evolution. Hence 
strong interaction-induced atom-number correlations can 
develop between the atomic wavepackets. 

Typical in-phase and out-of-phase bound states of two 
gap solitons in the middle of spectral gap are shown in 
Fig. 13 These states contain the same number of atoms 
for the same /i, which equals to exactly twice the num- 
ber of atoms in a single soliton. Remarkably, the internal 
noise correlation properties are the same for both pairs, 
which results in their identical optimal number squeezing 
ratio, see Fig. |SJa). Soliton entanglement can be quanti- 
fied by inter-soliton atom number correlations which can 
be calculated according to Eq. with An^.j being per- 
turbations of the atom number operator of two solitons, 
and indices {«, j} = {1,2} numbering individual solitons. 



The evolution of the correlation parameter C12 shows 
dramatic difference between the in- and out-of-phase soli- 
ton pairs. As follows from Fig. |SJb), the atom number 
correlation parameter for the in-phase soliton pair is neg- 
ative while it is positive for the out-of-phase pair. The 
difference in the entanglement behavior can be under- 
stood by examining the interaction energy of the solitons 
in the pair £'12 = / E\2(x)dx^ with the energy density: 

= E(^$^ + ^wv'.^. + mm] + 4^fV',). 

Interaction energy, Eyi depends on overlapping tail struc- 
ture of the stationary profiles ■(/'i,2(2;) corresponding to 
mean- field solitons ^'o. This quantity is small but posi- 
tive, = 3.1 X 10~^, for the out-of-phase soliton pair, 
and negative, = —3.4 x 10~^, for the in-phase pair. 
This results in positive and negative atom-number corre- 
lation parameter for the respective pairs. Small difference 
in the interaction energies is responsible for the slight 
asymmetry in correlations development seen in Fig. |S1 
(b). As a result, longer evolution times are required for 
the development of maximum atom number correlations 
between the solitons in an in-phase pair. 

In conclusion, we have investigated the effect of the 
periodic potential on quantum fluctuations and entan- 
glement properties of gap solitons in repulsive BEC con- 
fined by an optical lattices. We have revealed that the 
quadrature squeezing of gap solitons can be significantly 
enhanced compared to the matter-wave solitons in the 
lattice-free case. This enhancement is a strong effect that 
can potentially be detected in experiments on BEC gap 
solitons in one-dimensional optical lattices. We would 
like to emphasize that the basic results of our analysis 
can be useful in the study of the bandgap effects on the 
quantum squeezed states in other fields, such as quantum 
optics and gap solitons in photonic crystals. 
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